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Pyrolysis  and  Boundary  Layer  Combustion  of  a  Non- 
Charring  Solid  Plate  Under  Forced  Flow 

1.0  INTRODUCTION 

Burning  rate  is  a  measure  of  flammability  of  a  material  and  is  a  critical  parameter  in  fire 
threat  evaluation  and  for  designing  fire-fighting  systems.  Burning  vertical  walls  and  ceilings 
often  form  boundary  layer  flames  in  practice  in  which  burning  rates  play  a  fundamental  role  in 
determining  flame  size  and  flame  spread  [1].  To  gain  insights,  theoretical  and  laboratory  studies 
often  employ  well  controlled  forced  convection,  instead  of  the  natural  convection,  and  form  non¬ 
spreading  boundary  layer  flames  on  a  non-charring  polymer  plate,  which  has  well  defined 
physical  properties  (e.g.,  polymethyl  methacrylate,  PMMA)  [2].  In  the  boundary  layer 
configuration,  gas  phase  transport  and  heat  flux  from  the  flame  to  the  surface  decrease 
significantly  with  distance  along  the  fuel  surface  as  the  boundary  layer  thickness  increases.  The 
heat  feedback  drives  the  pyrolytic  reactions  and  in-depth  heat  and  mass  transfer  in  the  condensed 
phase.  Due  to  the  large  thermal  mass  of  the  condensed  phase,  transport  occurs  at  a  slow  rate 
involving  a  long  transient  unlike  the  transport  in  the  gas  phase  [3].  Furthermore,  the  heat  flux 
may  also  depend  on  the  shape  of  the  surface,  which  may  change  significantly  with  time  within 
the  pyrolysis  zone.  Therefore,  the  pyrolysis  or  burning  rates  vary  with  distance  along  the  surface 
as  well  as  the  bum  time  and  are  not  clearly  understood.  Also,  the  classical  boundary  layer  flame 
approximations  may  not  be  accurate  within  the  pyrolysis  zone,  which  may  not  extend  to  a  large 
enough  distance  from  the  leading  edge  due  to  the  long  transient  created  by  the  small  values  of  the 
surface  heat  flux. 

Emmons  [4]  obtained  similarity  solutions  ofNavier-Stokes  (NS)  equations  by  making 
classical  boundary  layer  approximations,  which  consist  of  infinite  rate  kinetics,  unit  Lewis  Le, 
Prandtl  Pr  and  Schmidt  Sc  numbers,  constant  surface  temperature  and  composition.  He  obtained 

an  analytical  expression  that  showed  that  the  local  burning  rate  or  mass  loss  rate  m  varies 
inversely  with  the  square-root  of  distance  from  the  leading  edge  Xj  and  directly  with  the  square- 

root  of  the  air  velocity  U.  Therefore,  the  dimensionless  regression  rate  or  Nusselt  number  Nu  is 
proportional  to  Re1/2,  where  Re  is  the  Reynolds  number.  Due  to  the  simplicity  of  the  closed  form 
equation,  the  Emmons  classical  solutions  have  been  used  widely  and  provide  a  fundamental  basis 
for  understanding  combustion  of  liquid  and  solid  fuels.  Sibulkin  et  al.  [5]  obtained  similar 
solutions  for  buoyancy-driven  flow  instead  of  forced  flow  considered  by  Emmons  [4],  but 
included  the  effects  of  finite  rate  combustion  kinetics.  Chen  and  Tien  [6]  obtained  numerical 
solutions  of  NS  equations  at  a  small  neighborhood  (few  mm  @  U=100  cm/sec)  near  the  leading 
edge  of  a  forced  convection  boundary  layer  flame  without  some  of  the  boundary  layer 
approximations.  However,  the  solutions  focused  on  flame  attachment  at  small  values  of  the 
Reynolds  number,  Re=Ux,/v  <10,  rather  than  the  burning  rate  distributions.  Their  solutions 

showed  the  effect  of  decreasing  Damkohler  number  Da  on  flame  transition  from  a  “closed”  to  an 
“open”  flame.  Mao,  Kodama,  and  Femandez-Pello  [7]  presented  NS  solutions  for  mixed 
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convection  effects  on  thermal,  concentration  and  velocity  fields  in  a  boundary  layer  flame  rather 
than  the  burning  rate  distributions.  Kodama,  Miyasaka  and  Femandez-Pello  [8]  calculated  the 
effects  of  air  flow  velocity  and  oxygen  mass  fraction  on  extinction  distance  for  the  mixed 
convection  boundary  layer  flame  modeled  by  Mao  et  al.  [7].  They  compared  the  predicted 
extinction  lengths  with  those  measured  from  diffusion  flames  formed  by  injecting  heptane 
through  a  porous  plate  burner.  In  this  work,  we  obtain  solutions  of  steady-state  NS  equations, 
which  are  similar  to  those  considered  by  Mao  et  al.  but  include  the  effects  of  composition- 
dependent  properties  and  surface-pyrolysis  kinetics,  which  is  imposed  by  developing  a  multi- 
variable  iterative  method.  Our  first  goal  is  to  predict  the  burning  rate  profiles  along  the  fuel 
surface  without  the  boundary  layer  approximations  and  make  quantitative  comparisons  with  the 
classical  Emmons  solutions. 

PMMA  forms  a  thin  layer  of  melt,  which  unzips  almost  entirely  to  the  monomer  (methyl 
methacrylate)  by  a  stepwise  mechanism  [9].  The  monomer  rapidly  vaporizes  and  escapes  from 
the  superheated  polymer  melt.  Therefore,  PMMA  pyrolysis  has  been  commonly  viewed  as  a 
surface  process  leading  to  quasi-steady-state  burning.  In  the  boundary  layer  configuration, 
however,  the  heat  flux  from  the  flame  to  the  surface  is  low  at  large  distances  from  the  leading 
edge.  At  the  low  surface  heat  flux,  the  thermal  boundary  layer  thickness  in  the  condensed  phase 
is  significant  and  the  pyrolysis  can  occur  in-depth.  Also,  the  heat-up  time  is  large  at  low  surface 
heat  flux  due  to  the  large  thermal  mass  of  the  condensed  phase.  Experiments  [  10,1 1]  on  PMMA 
plate  exposed  to  uniform,  pre-specified,  radiant  flux  clearly  showed  that  the  heat  transfer  into  the 
condensed  phase  occurs  slowly  due  to  the  large  thermal  mass  and  phase  change  phenomena.  In 
these  experiments,  which  were  performed  without  combustion,  the  rate  of  in-depth  pyrolysis 
increased  with  time  before  reaching  a  steady  state.  Furthermore,  the  time  to  obtain  a  steady  rate 
increased  significantly  as  the  radiant  heat  flux  decreased. 

Mekki  et  al.  [12]  performed  experiments  on  wind-aided  flame  spread  underneath  a 
PMMA  plate,  which  was  fixed  to  the  ceiling  of  a  rectangular  channel.  Mekki  et  al.  measured  the 
total  specie  production  rate  and  inferred  that  the  total  pyrolysis  rate  increased  with  time.  They 
also  found  that  the  total  specie  production  rate  increased  linearly  with  the  length  of  the  pyrolysis 
zone  contrary  to  the  square-root  relationship  predicted  by  the  Emmons  theory.  Agrawal  and 
Atreya  [3]  generalized  Emmons  solutions  by  combining  it  with  a  one-dimensional  theory  for  the 
transient  heat  up  and  pyrolysis  of  a  PMMA  plate.  Agrawal  and  Atreya  [3]  showed  qualitative 
agreement  between  their  theory  and  the  measurements  of  Mekki  et  al  [12]  for  the  total  production 
rate  of  species.  However,  very  few  studies  focused  directly  on  the  burning  rate  distributions  and 
the  role  of  transient  pyrolysis,  which  remains  unclear  for  a  forced-flow  boundary  layer  flame. 

Zhou  and  Femandez-Pello  [13]  measured  local  regression  rate  distributions  for  wind- 
aided  flame  spread  over  a  PMMA  plate  at  a  large  (>15  cm)  distance  from  the  leading  edge.  The 
flame  spread  experiments  were  conducted  both  without  and  with  externally  imposed  turbulence 
on  the  bulk  flow.  The  dimensionless  regression  rates  Nu  were  found  to  be  proportional  to  Re1/2 
for  laminar  flow  in  qualitative  agreement  with  Emmons  solutions  at  large  distances  (15  cm  <x  < 

70  cm)  from  the  leading  edge.  Their  results  also  seem  to  be  consistent  with  the  data  discussed  by 
Ahmed  and  Faeth  [14]  for  buoyancy-driven  flames  formed  over  a  vertical  plate.  Agrawal  [15] 
performed  experiments  using  an  apparatus  similar  to  the  one  used  by  Mekki  et  al.  [12]  for  flame 
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spread  underneath  a  PMMA  surface.  Agrawal  [15]  measured  the  regression  rates  for  distances  (0 
<x  <40  cm)  from  the  leading  edge.  At  large  values  of  x|5  the  heat  flux  arriving  at  the  surface  is 

small,  and  one  might  expect  a  significant  in-depth  pyrolysis  and  a  long  transient  for  heat-up  of 
the  cold  PMMA  plate.  Therefore,  Agrawal  [15]  attempted  to  account  for  the  unsteady  transport 
in  the  condensed  phase  using  the  theory  of  Agrawal  and  Atreya  [3]  and  reported  the  surface  heat 
flux  profiles  qw,  which  are  time  independent.  Near  the  leading  edge  of  the  plate,  the  qw  profile 
exhibited  a  peak.  The  flux  increased  within  the  first  3  cm  from  the  leading  edge  of  the  plate  and 
decreased  with  the  increased  axial  distance  xr  The  steady-state  surface  flux  was  found  to 

decrease  as  x^0  8  rather  than  x^0  5  predicted  by  Emmons  at  large  values  of  xp  Unlike  the  flame 

spread  experiments,  Krishnamurty  and  Williams  [2]  established  quickly  (5  sec)  a  forced-flow 
boundary  layer  flame  by  igniting  the  entire  surface  of  a  PMMA  plate  and  reported  regression 
rates  for  different  concentrations  of  oxygen  in  the  bulk  flow  of  the  oxidizer.  The  regression  rates 
were  measured  for  a  fixed,  short  bum  time  (MOO  sec)  at  three  positions  along  the  length  (15  cm) 
of  the  plate.  The  data  were  found  to  be  in  good  agreement  with  Emmons  solutions  for  high 
oxygen  concentration.  At  ambient  air,  however,  the  measured  regression  rates  were  found  to  be 
significantly  lower  than  the  predictions  of  the  classical  theory.  The  exact  reasons  for  these 
discrepancies  among  the  data  and  the  classical  theory  have  not  been  clear. 

Emmons  solutions  and  the  NS  solutions  presented  here  assume  that  the  fuel  surface 
remains  flat  as  it  bums  at  a  steady  rate.  However,  the  surface  heat  flux  and  the  local  regression 
rates  are  non-uniform  along  the  PMMA  plate.  Therefore,  the  surface  regresses  non-uniformly 
and  an  initially  flat  surface  of  a  real  material,  such  as  PMMA,  develops  curvature  as  the  time 
progresses.  Indeed,  experiments  [16,17]  performed  in  our  laboratory  have  recently  shown  that 
the  surface  curvature  affects  local  heat,  mass  and  momentum  transport,  which  in  turn  can  have 
significant  effect  on  the  local  burning  rates  for  PMMA.  Furthermore,  Hegab  et  al.  [18]  have 
recently  obtained  self-consistent  solutions  to  describe  burning  of  a  sandwich  propellant  exposed 
to  quiescent  air  and  showed  significant  moving  boundary  effects  on  the  burning  rates. 

Therefore,  our  second  objective  is  to  make  quantitative  comparisons  of  the  steady-state 
NS  solutions,  and  Emmons  solutions  with  our  experimental  measurements  of  regression  rate 
distributions  [16,17]  at  different  bum  times  (5  to  20  min)  to  study  the  transient  effects  and  the 
validity  of  the  approximations  (flat  surface  and  steady  surface-pyrolysis).  The  experiments  were 
performed  on  non-spreading,  forced-flow,  boundary  layer  flames,  which  were  formed  quickly  by 
exposing  the  entire  surface  of  a  PMMA  plate  (10  cm)  to  uniform  radiant  flux  rather  than  allow 
co-current  flame  spread.  Therefore,  every  point  on  the  surface  was  exposed  to  the  flame  for  a 
prespecified  length  of  time,  which  was  varied  systematically  in  the  experiments  until  steady-state 
is  approached. 

2.0  APPROACH 

Our  approach  includes  both  numerical  modeling  and  laboratory  experiments.  We  will 
focus  here  mainly  on  the  theoretical  aspects  of  the  problem  and  provide  only  a  brief  description 
of  the  experiments.  Further  details  of  the  experiments  can  be  found  in  a  separate  NRL  memo 
report,  NRL/MR/6 180-02-8630  [17]. 
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2.1  Theoretical 

We  consider  laminar  flow  of  air  past  a  flat  solid  surface,  which  consists  of  a  leading  inert 
plate  (1.7  cm),  followed  by  a  PMMA  plate  (6.7  cm),  followed  by  another  inert  plate  (8  cm)  as 
shown  in  Fig  1.  The  leading  plate  establishes  the  no-slip  condition  and  a  Blassius  momentum 


Figure  1  Boundary  Layer  Flame  Formed  Over  a  PMMA  Plate 


boundary  layer.  The  lengths  of  the  inert  plates  were  chosen  to  minimize  the  effects  of  in-flow 
and  out-flow  boundary  conditions  on  the  combustion  region.  The  width  of  the  plate  is  assumed 
to  be  infinite.  The  air  enters  the  domain  at  a  uniform  velocity  parallel  to  the  plate  and  forms  a 
boundary  layer  over  the  surface.  The  nominal  Reynolds  number,  Re=Ux|/v  is  7000  (laminar)  or 
less,  where  Xj  is  the  distance  from  the  leading  edge  of  the  fuel  plate,  x=x-1.7,  where  x  is 
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measured  from  the  leading  edge  of  the  inert  plate.  The  nominal  value  of  Re  is  evaluated  using 
the  ambient  value  of  the  kinematic  viscosity,  v=0.157  cm2/sec.  Mao  et  al.  [19]  showed  that  the 
effects  of  buoyancy  and  forced  flow  on  combustion  scale  as  Gr/Re2  5  for  a  horizontal  fuel  surface 
rather  than  the  commonly  used  relationship  Gr/Re2  for  a  vertical  surface.  A  small  difference  in 
the  exponent  of  Re  has  a  large  effect  on  the  ratio  and  on  the  significance  of  buoyancy  effects. 
Therefore,  for  the  small  sample  length  (<10  cm)  and  for  U>60  cm/sec,  Gr/Re2'5  <0.18  and  the 
forced  convection  are  assumed  to  dominate  the  gravity  effects. 

Initially,  the  air  and  the  surface  are  kept  at  room  temperature,  and  the  entire  fuel  surface  is 
ignited  by  adding  external  energy  to  the  gases  above  the  plate  for  a  short  period  of  time.  Upon 
ignition,  the  PMMA  undergoes  pyrolysis  to  form  monomer  vapor,  which  mixes  with  the  air  by 
convection  and  diffusion.  The  fuel-air  mixture  undergoes  combustion  and  a  boundary  layer 
flame  is  established  above  the  fuel  plate.  The  unsteady,  full  NS  equations  describe  the 
development  of  the  flame  from  ignition  to  steady  state.  The  coordinates  x  and  y  are  fixed  to  the 
leading  edge  of  the  leading  inert  plate  and  represent  distances  along  the  surface  and 
perpendicular  to  the  surface,  respectively.  The  conservation  equations  are  given  in  Appendix  A. 

Combustion  of  the  methyl  methacrylate  monomer  is  assumed  to  take  place  in  a  single 
step  to  produce  stochiometric  amounts  of  carbon  dioxide  and  water  according  to 

C5H,02  +  602  =  5 C02  +  4 H20  .  (1) 

Therefore,  the  production  rates  of  CO2  and  H2O  and  consumption  rates  of  O2  are  proportional  to 
the  rate  of  combustion  of  fuel  Wk,  which  is  assumed  to  follow  second  order  reaction  kinetics 
given  by  Seshadri  and  Williams  [20]. 

Wk=kcCk  C02,  (2) 

where  kc  is  the  kinetic  constant  in  cm3/gmole-sec.  It  is  given  by 

kc  =  2.3xl016  exp(-43000  /  RT ) ,  (3) 

where  T  and  R  are  the  gas  temperature  and  the  gas  constant,  respectively.  The  energy  release  rate 
per  unit  volume  HRRV  is  given  by 

HRRv  =  AHcWk  /  Av ,  (4) 

where  AHC  (2530  KJ/gmole)  and  Av  are  the  heat  of  combustion  and  Avogadro  number, 
respectively. 

Binary  diffusion  coefficients  Dki,  viscosity  |iki  and  conductivity  A.ki  are  calculated  for  each 
specie  from  power  law  correlation  of  Chapman-Enskog  equations  as  functions  of  temperature  as 
shown  in  Appendix  A  The  mole  fraction  of  fuel  vapor  in  the  combustion  gases  is  typically  much 
smaller  than  1 .0.  The  properties  of  n-pentane  gas  are  used  for  the  pure  monomer  vapor  after 
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correcting  for  the  molecular  weight  difference,  which  is  about  30%.  Specific  heats  of  pure 
species  are  taken  from  standard  tables.  The  properties  of  gas  mixture  are  evaluated  from  the  pure 
specie  values  using  the  mixture  rules  described  in  Appendix  A 

As  the  heat  from  the  flame  is  transferred  to  the  polymer  surface,  the  polymer  first  softens 
and  begins  to  melt  around  100  °C  [9].  Upon  further  heating,  the  melt  begins  to  pyrolyze  around 
250  °C  and  unzips  to  form  the  monomer.  Therefore,  the  monomer,  which  has  a  normal  boiling 
point  of  about  101  °C,  instantly  evaporates  from  the  super-heated  polymer  melt,  and  forms 
bubbles,  which  are  ejected  to  the  surface.  As  a  result  the  melt  is  mostly  composed  of  the  polymer 
and  is  very  viscous.  The  thickness  of  the  melt  layer  depends  on  the  heat  flux  to  the  surface.  At 
high  surface  heat  flux,  which  occurs  near  the  flame  attachment  point,  the  melt  layer  is  barely 
visible  due  to  high  temperature  gradients,  and  the  pyrolysis  reaction  may  be  assumed  to  occur  at 
the  surface.  However,  as  x{  increases  and  qw  decreases,  the  melt  layer  becomes  thick  due  to 

small  temperature  gradients,  and  pyrolysis  may  occur  in  depth  at  an  unsteady  rate.  However,  due 
to  the  small  length  of  the  PMMA  plate  considered  in  this  analysis,  we  assume  that  the  melt  layer 
is  very  thin  and  the  pyrolysis  reaction  occurs  essentially  at  the  solid  surface  at  a  steady  rate 
following  an  ablative  step; 

[{-C.Hfi-X  L„  n[C5Hs02U,  (5) 

The  rate  of  pyrolysis  is  determined  both  by  the  rate  of  heat  transfer  to  the  polymer  and  the 
intrinsic  kinetics  of  the  pyrolysis  reaction.  Therefore,  in  the  limit  of  infinite  heat  transfer,  the 
pyrolysis  rate  is  determined  by  the  intrinsic  chemical  kinetics  of  the  pyrolysis,  which  are 
measured  under  isothermal  conditions  to  eliminate  resistance  to  heat  transport.  On  the  other 
hand,  in  the  limit  of  infinite  rate  chemical  kinetics,  the  pyrolysis  rate  is  determined  by  the  rate  of 
heat  transfer  to  the  polymer.  We  assume  that  the  intrinsic  chemical  kinetics  follows  a  zeroth 
order  kinetic  expression  given  by  Arisawa  and  Brill  [21].  Arisawa  and  Brill  [21]  performed 
experiments  with  very  thin  films  of  anionically  polymerized,  high-molecular- weight  (131,000) 
PMMA  at  different  values  of  the  preset  temperature  Ts  (653-733  K).  The  temperature  was  kept 
uniform  (+/-  2  °K)  across  the  film  in  their  experiments  by  supplying  heat  to  compensate  for  the 
cooling  due  to  endothermic  pyrolysis.  Therefore,  they  provide  intrinsic  measurements  of  the 
kinetics  that  depend  only  on  the  temperature,  which  was  calibrated  against  compounds  with 
standard  melting  points.  The  temperature  is  kept  uniform  across  the  film  in  their  experiments 
and,  therefore,  provide  intrinsic  measurements  of  the  kinetics  that  depend  only  on  the 
temperature  and  independent  of  heat  transfer  limitations.  Their  rate  constants,  which  are  in  1/sec, 
were  multiplied  by  the  thickness  (83  pm)  of  the  film  and  the  polymer  density  (ps=l  .19  gm/cm  ) 

to  obtain  an  expression  for  the  surface  mass  flux  or  the  pyrolysis  rate,  m  in  gm/cm2sec,  which  is 
given  by 

r0ps  =m  =  8x1 017  exp(-66000/  RTs )  .  (6) 
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Equation  (6)  shows  an  exponential  relationship  between  the  pyrolysis  rate  and  the  surface 
temperature,  which  varies  along  the  surface.  Therefore,  an  order  of  magnitude  change  in  the 
pyrolysis  rate  results  in  a  small  change  of  35  K  in  surface  temperature  due  to  the  high  activation 
energy  (66  Kcal/gmole  or  2.76  xlO5  KJ/Kgmole).  Here,  rQ  is  the  local  regression  rate  of  the 

surface.  Significantly  lower  values  than  66  Kcal/gmol  for  Ea  were  reported  by  others  [2,10]. 
Kashiwagi  et  al.  [22]  showed  that  the  activation  energy  is  generally  much  larger  (56  Kcal/gmole) 
in  the  absence  of  oxygen  at  the  PMMA  surface  than  that  (15  Kcal/gmole)  for  oxidative 
degradation.  Our  simulations  show  very  little  oxygen  concentration  on  the  PMMA  surface  and 
also  rQ  to  be  less  sensitive  to  Ea  as  discussed  later. 

At  the  PMMA  surface,  the  pyrolysis  rate  must  equal  the  sum  of  convective  and  diffusive 
fluxes  for  the  monomer  vapor  and  is  stated  by 


m  =  m Xt(t,x, 0)  -  D^pdX^t^fi)  /  dy.  (7) 

The  net  flux  for  oxygen,  nitrogen,  carbon  dioxide  and  water  vapor  at  the  surface  is  zero  and  is 
given  by 

m  Xk  (t,  x,0)  -  DkmpdXk  (t,  *,0)  /  dy  =  0,  (8) 

where  is  the  mole  fraction  of  the  specie  k.  As  the  boundary  layer  thickness  increases  along 

the  plate,  both  the  diffusive  flux  and  m  decrease.  Therefore,  equation  (8)  requires  that  Xk  vary 
along  the  surface.  The  heat  transfer  from  the  flame  to  the  surface  should  equal  the  heat  needed 
to  pyrolyze  the  solid  and  is  given  by 


mQP=  KdT (6  *>0)  I  dy.  (9) 


Here,  Qp  ( 1.62xl03  KJ/Kg)  is  the  endothermic  heat  of  pyrolysis  for  PMMA.  In  writing  equation 
(9),  net  radiative  heat  feedback  is  neglected  compared  to  the  gas  phase  heat  transfer  by 
conduction  and  convection  for  the  small  length  of  the  PMMA  plate  considered  in  this  analysis. 
The  net  radiation  is  the  result  of  radiative  heat  feedback  to  the  surface  and  re-radiative  cooling 
from  the  surface,  which  may  be  comparable  due  to  high  emissivity  of  black  PMMA  surface  [2]. 
The  in-depth  heat  transfer  and  pyrolysis  are  expected  to  become  significant  at  increasing 
distances  from  the  leading  edge  of  the  plate  as  discussed  later,  and  equation  (9)  becomes  less 
accurate.  In-depth  heat  transfer  is  also  assumed  to  be  small  for  the  small  length  of  the  plate 
considered  in  this  study.  The  velocity  components  at  the  PMMA  surface  are  given  by 

jfl 

K  =  ~T+ — —  ,u(t,x,0)  =  0  .  (10a, b) 

p{t,x,  0) 
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The  dimensionless  regression  rate  Nu  is  defined  as 


where  the  thermal  conductivity  of  gas  at  the  PMMA  surface  Xw=0. 00052 W/m2K  is  taken  as  the 
value  for  air  at  700  °K.  Also,  the  temperature  drop  across  the  flame  AT=1200  °K  was  used  since 
it  did  not  change  significantly  with  x  for  the  relatively  small  length  (<10  cm)  of  the  fuel  plate 
considered  in  this  study. 

At  low  air  velocities,  the  flame  attaches  close  to  the  leading  edge  of  the  PMMA  plate. 
Therefore,  the  heat  transfer  characteristics  of  the  leading  plate  in  front  of  the  sample  could  play  a 
role  in  the  attachment  dynamics.  Therefore,  heat  loss  through  the  leading  plate  is  assumed  to 
follow  Newton’s  law  of  cooling  with  a  characteristic  heat  transfer  coefficient  (  h=173  W/m2K). 
No-slip  condition  was  applied  along  the  entire  surface.  The  no-slip  was  applied  only  on  a  small 
portion  of  the  leading  plate  to  accurately  represent  the  experimental  conditions.  The  post- flame 
plate  temperature  is  kept  uniform  at  the  trailing  edge  temperature  of  the  PMMA  plate.  At  the  top 
of  the  computational  domain,  air  is  allowed  to  flow  axially  with  slip.  The  boundary  conditions 
for  the  inert  plates,  inlet  and  outlet  are  also  given  in  Appendix  A. 

Due  to  high  thermal  mass  of  the  condensed  phases,  the  pyrolysis  and  transport 
phenomena  occur  at  a  very  slow  rate,  unlike  the  phenomena  in  the  gas  phase.  Therefore,  the 
surface  is  assumed  to  remain  flat  despite  non-uniform  heat  feedback  and  regression  of  the 
surface.  The  flat  plate  assumption,  together  with  the  surface-pyrolysis  assumption  leads  to  quasi¬ 
steady-state  in  which  all  of  the  dependent  variables  in  the  conservation  equations  and  the 
burning  rate  distributions  become  independent  of  time  at  a  short  time  (<  2  sec)  after  the  ignition. 
Therefore,  equations  (6-10)  are  applied  at  y=0.  The  implications  of  these  assumptions  on  the 
comparison  between  the  steady-state  NS  solutions  and  the  experiments  will  be  discussed  later. 


2.2  Numerical 

There  are  several  unknowns,  i.e.,  p,  u,  v,  P,  E,  and  Ck  in  the  gas  phase,  as  well  as 
temperature,  specie  concentrations,  and  vertical  velocity  distributions  along  the  surface  of  the 
PMMA  plate.  All  of  these  are  calculated  as  functions  of  position  x  and  y  and  time  t.  One 
typically  specifies  inlet  velocity  uo,  density  po,  composition  Cko,  and  temperature  To  (or  energy 
Eo)  of  the  air.  The  equations  are  discretized  by  a  finite  volume  method  on  an  Eularian  domain 
containing  about  192x144  finite  volume  cells  in  x  and  y  directions,  respectively.  The  cells  are 
clustered  closely  near  the  leading  edge  of  the  flame,  which  is  not  known  a  priori.  At  large 
velocities,  the  flame  attaches  downstream  and  the  cell  placement  is  corrected  accordingly.  The 
minimum  cell  size  inside  the  flame  is  0.2  mm  x  0.2  mm.  Therefore,  at  the  very  leading  edge,  as 
Xj  becomes  comparable  to  0.2  mm,  the  solutions  presented  here  become  less  accurate.  The  cells 

are  stretched  in  both  directions  from  the  leading  edge.  The  convective,  diffusive,  combustion 
and  radiation  contributions  are  solved  by  separate  algorithms  using  separate  time  steps,  which  are 
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typically  either  smaller  or  equal  to  the  global  time  step.  Diffusion  and  reaction  processes  are 
subcycled  within  a  global  time  step  and  cumulative  contributions  from  each  process  are 
calculated.  These  terms  are  then  coupled  together  by  time-step  splitting  as  discussed  by  Oran 
and  Boris  [23].  The  global  time  step  is  chosen  such  that  the  diffusion  processes  are  not 
subcycled  by  more  than  three  times.  The  global  time  step  is  typically  of  the  order  of  10  and  50 
msec  during  and  after  ignition,  respectively. 

Ignition  along  the  entire  length  of  the  plate  was  achieved  by  adding  external  energy  to 
several  number  of  rows  of  cells  above  the  plate  for  about  msec.  During  the  ignition  time,  the 
PMMA  surface  temperature  was  prespecified  at  700  K  to  ensure  that  the  fuel  vapor  is  ejected 
continuously  from  the  surface  forming  combustible  mixture  above  the  plate.  Therefore,  it 
simulates  qualitatively  the  radiative  ignition  employed  in  the  experiments,  which  are  discussed 
later.  The  magnitudes  of  ignition  time  scales  are  much  smaller  in  the  numerical  simulations  than 
those  employed  in  the  experiments.  Therefore,  the  transient  behavior  is  expected  to  differ 
between  the  simulations  and  the  experiments  for  the  gas  phase.  However,  the  ignition  method 
has  little  effect  on  the  steady  state  results  for  the  gas  phase,  which  is  the  main  focus  of  the 
computations. 

Convection  equations  are  solved  by  the  high-order  implicit  algorithm,  Barely  Implicit 
Correction  to  Flux-Corrected  Transport  (BIC-FCT),  which  was  developed  by  Patnaik  et  al.[24]  to 
solve  the  convection  part  of  the  equations  for  slow  flows.  The  Flux-Corrected  Transport  (FCT) 
is  an  explicit  scheme  for  the  solution  of  hyperbolic  equations.  Therefore,  it  requires  very  small 
time  steps  dictated  by  sound  speed  following  Courant-Friedrichs-Lewy  condition.  BIC-FCT  was 
developed  based  on  the  suggestion  by  Casuli  and  Greenspan  [25]  that  only  the  pressure  in 
momentum  equation  and  velocity  in  energy  equations  needs  to  be  treated  implicitly  to  increase 
the  time  step  to  that  based  on  fluid  velocity  rather  than  the  sound  speed.  This  method  involves 
solution  of  an  elliptic  equation  for  pressure  which  is  solved  using  MGRID  algorithm  discussed  in 
detail  by  DeVore  [26].  Typically,  the  courant  number  based  on  maximum  fluid  velocity  and 
minimum  cell  size  is  smaller  than  0.3. 

Diffusion  equations  are  discretized  using  explicit  finite  difference  approximations.  The 
time  steps  for  the  subcycle  are  calculated  based  on  minimum  cell  size  and  maximum  diffusivity 
in  the  domain.  If  these  time  steps  are  more  than  one  third  of  the  global  time  step,  then  the  global 
time  step  is  reduced.  At  the  end  of  each  subcycle,  diffusion  contributions  to  specie 
concentrations,  temperature,  and  momentum  are  updated  for  mass,  heat  and  momentum, 
respectively. 

Reaction  contribution  is  described  by  a  single  nonlinear  stiff  ordinary  differential 
equation  and  is  solved  by  using  Burlich-Stoer  algorithm  with  Rhomberg  interpolation  as 
described  by  Press  et  al.  [27].  Unlike  the  diffusion  processes,  the  reaction  time  steps  are 
evaluated  based  on  local  values  of  the  variables  rather  than  global  maximum/minimum  values  in 
the  domain.  The  time  steps  are  calculated  cell  by  cell  from  the  local  values  of  rate  constant, 
specie  concentrations  and  temperature  and  vary  widely  across  the  domain.  Typically,  the  time 
steps  can  be  hundreds  of  times  smaller  than  the  global  time  step,  especially  during  ignition  near 
the  leading  edge  of  the  porous  plate.  After  ignition,  the  time  step  increases  by  a  factor  of  3-10  in 
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the  hottest  regions.  Therefore,  the  reaction  calculations  are  subcycled  to  a  varying  degree  within 
the  domain.  Specie  concentrations  are  updated  at  the  end  of  each  subcycle.  A  very  significant 
cost  of  computation  is  incurred  in  reaction  calculations. 

Calculation  of  the  burning  rates  require  special  attention.  Unlike  the  gas  phase,  which  is 
treated  as  a  continuum,  the  PMMA  surface  is  treated  as  plane  of  discontinuity  on  which 
equations  (6-10)  are  imposed  by  a  multi-variable  fixed  point  iterative  method.  In  this  method,  at 
each  time  step,  new  values  of  rh  ,  Xk  and  Ts  are  evaluated  from  equations  (6)  to  (10),  in  which 
old  values  are  used  to  calculate  the  derivatives.  Solution  of  the  pyrolysis  equation  (6)  needs 
special  attention.  An  implicit  form  of  the  pyrolysis  equation  was  derived  by  substituting 
equation  (6)  into  (9)  and  was  solved  by  Newton-Raphson  iteration  with  bisection  to  obtain  the 
new  value  of  Ts.  Therefore,  within  each  fixed-point  iteration,  only  Ts  was  solved  by  the  Newton- 
Raphson  method  (inner  iteration),  with  both  rh  and  Xk  fixed.  If  the  new  values  of  rh  ,  Xk  and  Ts 
are  not  within  a  prespecified  tolerance  of  the  values  from  previous  fixed  point  (outer)  iteration, 
the  Navier-Stokes  equations  are  re-solved  with  the  new  values  for  the  boundary  conditions  at  the 
surface.  Thus,  rh  ,  Xk  and  Ts  are  iterated  until  convergence  is  achieved  for  every  cell  along  the 
polymer  surface,  before  proceeding  to  the  next  time  step.  The  fixed-point  iterations  are  also  used 
to  implement  the  mixed  boundary  condition,  equation  (A34),  for  calculating  the  surface 
temperature  of  the  leading  plate,  without  the  Newton-Raphson  inner  iterations.  A  similar  method 
was  implemented  by  Ananth  et  al.  [28,29]  for  studying  a  boundary  layer  flame  formed  over  a 
porous  plate  without  the  inner  Newton-Raphson  iterations. 

2.3  Emmons  Boundary  Layer  Theory 

One  may  evaluate  the  accuracy  of  the  BL  approximations  by  comparing  the  NS  solutions 
with  the  classical  solutions.  Emmons  recast  the  steady-state  form  of  the  NS  equations  given  in 
Appendix  A,  equations  (6-10),  and  derived  an  analytical  solution  by  using  classical  boundary 
layer  approximations  described  by  Williams  [30],  For  large  Re,  which  occur  at  high  velocities  or 
large  values  of  x|}  axial  diffusion  terms  in  NS  equations  become  negligible  compared  to  lateral 
diffusion.  Also,  assuming  Le=Pr=Sc=land  the  combustion  kinetics  occur  at  infinite  rate  (flame 
sheet  approximation),  the  NS  equations  were  reduced  to  a  homogeneous  parabolic  equation, 
which  was  solved  for  the  burning  rate  using  similarity  method.  In  this  method,  intrinsic  kinetics 
of  surface  pyrolysis  was  assumed  to  occur  at  infinite  rate  so  that  the  pyrolysis  rate  is  determined 
by  the  heat  feedback.  Equation  (6)  was  replaced  by  constant  surface  temperature  boundary 
condition,  which  was  prespecified.  The  boundary  layer  solution  is  given  by 


f 


roPs  ~  PJW  ~ 


V 


2xl 


\ 


J 


1/2 


(11) 


Here,  subscript  w  and  oo  represent  values  at  the  the  solid  surface  and  in  the  bulk,  respectively. 
The  values  of  f(0)  are  represented  in  analytical  form  by  Glassman  [31]  as  shown  below 
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(12) 


ln(l  +  B) 

2.6B°U 

where  B  is  the  Spaulding’s  mass  transfer  number.  Equation  (1 1)  can  be  written  in  dimensionless 
form  as 


X  X  AT 

w  w 


PvKcQp 


(-/(0))VRe. 


(13) 


For  PMMA,  the  Spalding’s  mass  transfer  number  B  is  given  by  Kanury  [32]  and  Chen  and  Tien 
[6]  as  1.3,  which  is  substituted  into  equation  (12)  to  get  f(0)=  -  0.3.  The  value  of  f(0)  is  slightly 
smaller  than  the  value  (-0.397)  given  by  Krishnamurthy  and  Williams  [2],  Substituting  equation 
(12)  into  (13),  one  obtains  for  the  dimensionless  regression  rate  Nu  for  PMMA  (B=1.3) 


(14) 


An  expression  for  the  flame-standoff  distance  8  may  be  obtained  from  equation  (14)  by 
making  linear  approximation  for  the  temperature  gradient  at  the  surface  to  obtain 


i  _K  l 

Nu  hx )  x,  V  Re 


(15) 


Equation  (14)  shows  clearly  that  Nu  depends  only  on  a  single  parameter  Re.  It  increases 
with  square-root  of  U  and  decreases  with  square-root  of  x.  Therefore,  as  x  approaches  zero  as 
one  approaches  the  leading  edge,  Nu  reaches  a  very  high  value  due  to  the  neglect  of  axial 
diffusion  terms  in  the  NS  equations.  Indeed,  we  will  show  that  the  axial  diffusion  terms  remain 
significant  even  at  significant  distances  from  the  leading  edge.  Furthermore,  the  infinite  rate 
kinetics  approximation  in  the  NS  equations  leads  to  flame  “attachment”  at  the  leading  edge  of  the 
PMMA  plate  independent  of  U.  Full  NS  equations,  however,  predict  flame  “attachment”  at  some 
distance  from  the  leading  edge,  which  is  affected  by  the  competition  between  convection  and  the 
finite  rate  combustion  kinetics  as  discussed  later. 

2.4  Experimental 

Experiments  were  performed  on  non-spreading  boundary  layer  flame  formed  over  a 
relatively  short  length  of  black  PMMA  sample  (7.5  cm  x  9.5  cm  x  2.5  cm).  The  combustion 
experiments  were  conducted  for  different  values  of  the  bum  time  to  record  the  transient  effects. 
Experimental  details  are  given  elsewhere  [16,17,33-35]  and  only  a  brief  description  is  given  here 
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for  completeness.  The  geometry  is  similar  to  that  shown  in  Figure  1.  A  high  molecular  weight, 
black  PMMA  sample  was  placed  in  a  cup  holder,  which  has  extension  plates  (lips)  to  minimize 
the  edge  effects.  Thin  quartz  plates  are  placed  all  around  the  fuel  plate  to  minimize  excess 
burning  of  the  edges.  The  sample  is  placed  at  the  exit  of  a  channel  rather  than  inside  the  channel 
to  facilitate  temperature  mapping  of  the  flame.  Air  was  forced  out  of  the  channel  at  a 
prespecified  center  velocity.  The  velocity  distributions  across  the  air  jet  was  found  to  be  flat  and 
did  not  change  significantly  along  the  length  of  the  sample.  The  entire  PMMA  surface  was 
exposed  to  a  large  radiant  panel  typically  for  less  than  30  sec.  The  sample  surface  is  melted  and 
pyrolyzed  as  evidenced  by  the  formation  and  ejection  of  bubbles.  The  fuel  vapors  are  ignited  to 
quickly  form  a  boundary  layer  flame.  During  combustion,  a  thin  film  of  melt  is  formed  at  the 
surface  with  intense  bubbling  near  the  leading  edge  and  relatively  slow  bubbling  far  from  the 
leading  edge.  The  bubbles  were  very  small  and  barely  visible  near  the  leading  edge  and  increase 
in  size  with  xr  The  width  of  the  flame  is  much  larger  than  the  boundary  layer  thickness  so  that 

the  flame  can  be  approximated  as  being  two  dimensional. 

After  about  1  min  from  the  time  of  ignition,  five  fine  thermocouples  are  brought  into  the 
flame  to  map  the  temperature  profiles  across  the  flame  at  five  different  locations  Xj  along  the 

plate.  The  thermocouple  readings  are  reported  after  a  radiation  correction.  The  temperatures 
were  recorded  at  different  distances,  Y,  measured  from  the  initial  position  of  the  surface  rather 
than  the  regressing  surface.  The  sample  was  allowed  to  bum  for  a  prespecified  amount  of  time 
before  the  flame  was  extinguished,  and  the  sample  was  allowed  to  cool  down.  The  change  in  the 
thickness  of  the  plate  due  to  burning  was  measured  with  a  micrometer  at  different  distances  from 
the  leading  edge  xr  Regression  rate  is  the  change  in  the  thickness  of  the  sample  divided  by  the 

bum  time.  The  boundary  layer  formed  was  smooth  except  near  the  sides  and  near  the  trailing 
edge  of  the  sample  due  to  buoyancy  effects.  Measurements,  however,  were  made  well  within  the 
smooth  region  and  along  the  center  line  of  the  sample.  Some  of  the  experimental  results  are 
summarized  below  [16,17,33-35]: 

(1)  At  the  end  of  ignition,  the  PMMA  surface  was  flat.  As  the  time  progresses,  a  step  was 
formed  near  the  flame  attachment  point.  The  step  size  increased  in  depth  and  width  with  t.  In 
the  valley,  the  surface  to  flame  distance  increased  with  time  as  indicated  by  a  single 
thermocouple  measurements. 

(2)  At  a  fixed  value  of  the  bum  time,  the  regression  rate  rQ  increased  with  xt  and  reached  a 
maximum  near  the  flame  “attachment  point”  followed  by  a  decrease  in  rQ  with  xr 

(3)  rQ  decreased  significantly  with  bum  time  t  within  the  valley  at  short  distances  from  the 
leading  edge  xr  On  the  other  hand,  rQ  increased  with  t  at  large  values  of  and  was  often 

associated  with  increased  activity  of  bubbling  near  the  trailing  section  of  the  sample  plate.  As  U 
is  increased,  bubbling  activity  in  the  trailing  section  increased  significantly. 

(4)  At  large  values  of  U,  the  flame  attached  downstream  of  the  leading  edge  of  the  fuel  plate 
where  the  valley  was  formed.  As  the  time  progressed,  the  flame  spread  very  slowly  toward  the 
leading  edge  of  the  plate  while  the  valley  got  wider  in  size. 
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(5)  The  surface  was  semi-solid  and  the  melt  layer  was  barely  noticeable  near  the  leading  edge. 
The  melt  layer  thickness  appeared  to  increase  with  xr  The  melt  temperature  decreased  by  ~  45  K 
along  the  length  of  the  sample  as  indicated  by  pressing  a  thermocouple  onto  the  surface. 


3.0  RESULTS  AND  DISCUSSION 

Once  the  geometry,  the  combustion  kinetics  and  the  pyrolysis  kinetics  are 
specified,  one  only  needs  to  specify  the  air  velocity  to  obtain  solutions  of  the  governing 
equations.  Solutions  describe  the  temperature,  specie  mole  fraction,  velocity  and  density  profiles 
as  functions  of  x,  y  and  t  [36-41].  Also,  the  values  of  rh  ,  Xk  and  Ts  along  the  surface  are 
predicted.  Simulations  were  performed  until  steady  state,  which  is  achieved  in  less  than  a  second 
at  the  air  velocities  employed  in  our  experiments.  However,  our  goal  is  to  discuss  the  steady- 
state  values  for  gas-phase  temperature  profiles  and  burning  rates,  the  quantities  that  were 
measured  in  the  experiments  for  relatively  short  length  (6.5  cm)  of  the  PMMA  plate. 

Figure  2a  shows  typical  temperature  contours  for  inlet  air  velocity,  U,  of  60  cm/sec.  The 
combustion  reactions  occur  mostly  in  the  high  temperature  region,  which  is  represented  by  one  of 
the  dark  bands,  and  is  referred  to  as  the  “flame”  or  “combustion  zone”.  Near  the  leading  edge, 
fresh  fuel  comes  in  contact  with  fresh  air  and  the  combustion  rates  are  very  high.  Downstream, 
the  products  dilute  the  gases  and  the  combustion  rates  decrease  sharply  within  a  short  distance 
from  the  leading  edge  as  indicated  by  the  heat  release  rate  contour  (HRRV)  shown  as  a  dark  line 
labeled  “DEDtrxn”  for  lxlO9  ergs/cm3sec  in  Figure  2a.  The  rates  of  heat  and  mass  transfer 
between  the  flame  and  the  surface  are  also  large  near  the  leading  edge  due  to  steep  gradients,  and 
they  decrease  with  distance  from  the  leading  edge  as  the  boundary  layer  thickness  increases. 
Therefore,  near  the  leading  edge,  the  transport  and  combustion  rates  are  expected  to  be  closely 
coupled.  The  competition  between  convective  transport  and  reaction  kinetics  is  represented  by 
Damkohler  number  Da=  reaction  rate/convection  rate=  kpxi/U,  where  kp  =  kcCo2  is  the  psuedo- 
first-order  rate  constant.  Here,  C02  is  taken  as  the  ambient  concentration  of  oxygen  in  the  bulk. 

At  small  air  velocities,  Da  is  high  and  the  gases  get  heated  up  within  a  short  distance  from  the 
leading  edge  x,  by  the  high  rates  of  heat  release.  Therefore,  the  flame  “attaches”  near  the  leading 

edge  of  the  plate.  As  U  is  increased,  Da  decreases,  and  the  flame  “attaches”  at  increasing  value 
of  Xj  as  shown  in  Figure  2b.  Figures  2a  and  2b  also  show  that  the  flame  temperature  in  the 

combustion  zone  has  not  decreased  significantly  along  the  length  of  the  PMMA  plate.  Clearly, 
as  U  increses,  the  heat  release  rate  contour  (DEDtrxn)  increases  in  length  indicating  the  increased 
length  of  the  flame. 
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3.1  NS  for  Distributions  Along  the  Solid  Surface 

Figure  3  shows  the  NS  solutions  for  the  surface  temperature  profiles  along  the  length  of  the 
sample  at  different  air  velocities.  As  one  approaches  the  flame  “attachment”  points,  shown  in 
Figures  2a  and  2b,  the  surface  temperature  raises  sharply  to  a  maximum  and  then  decreases  with 
Xj  as  the  boundary  layer  gets  thicker  reducing  the  heat  feedback.  The  surface  temperature  changes 

by  about  40  K  over  6.5  cm  length  of  the  plate.  This  is  qualitatively  consistent  with  our 
experiments  on  PMMA  under  boundary  layer  combustion.  During  combustion,  the  surface  of 
PMMA  appears  to  be  covered  with  a  froth  layer  consisting  of  polymer  melt  and  monomer 
bubbles  formed  by  the  unzipping  of  the  polymer.  It  was  observed  that  the  bubbling  was  intense 
and  bubble  size  was  small  near  the  leading  edge.  At  large  values  of  x,,  the  bubbling  was  slow 

and  bubble  size  was  large.  The  rate  of  bubble  formation  should  increase  with  increased  rate  of 
pyrolysis,  which  occurs  at  high  melt  temperature.  At  high  melt  temperature,  the  superheat  is 
higher  and  the  bubble  size  may  be  expected  to  be  small.  The  experiments  seem  to  show  that  both 
the  bubble  size  increases  and  the  bubbling  rate  decreases  with  xr  Since  the  surface  is  not  a  well- 

defined  interface  between  solid  and  gas  phases,  it  is  difficult  to  measure  precisely  the 
temperature  changes  along  the  surface.  However,  we  attempted  to  measure  the  surface 
temperature  by  dipping  a  fine  thermocouple  into  the  froth  layer  on  the  surface  at  different 
distances  from  the  leading  edge  xf  We  estimate  the  temperature  drop  from  x  =0  to  x=9  cm  is 

about  40  K,  which  is  comparable  to  those  measured  in  our  experiments  and  by  Vovelle  et  al.  [10] 
in  their  experiments  on  flameless  pyrolysis  of  PMMA.  Vovelle  et  al.[10]  measured  a  change  of 
30  K  as  the  heat  flux  changed  from  3.0  W/cm2  to  1.5  W/cm2.  However,  the  values  of  Ts  and  its 
distribution  are  determined  mainly  by  the  kinetic  constants  in  the  pyrolysis  kinetics  expression 
given  by  equation  (6).  A  variation  of  these  constants  seem  to  affect  mainly  Ts  rather  than  the 
burning  rate  distributions.  Figure  3  also  shows  that  the  increased  U  moves  the  flame 
“attachment”  point  downstream  and  the  maximum  in  Ts  occurs  downstream.  Furthermore, 
increase  in  U  increases  Ts  slightly  at  large  values  of  xr 

The  specie  concentration  profiles  along  the  PMMA  surface  are  shown  in  Figures  4-6.  As  Ts 
decreases  along  the  surface,  the  rate  of  formation  of  MMA  decreases  as  indicated  by  the 
decrease  in  its  concentration  shown  in  Figure  4.  The  mole  fraction  of  MMA  in  Figure  4  is  close 
to  the  value  of  the  uniform  surface  concentration  used  in  the  calculations  of  Mao  et  al.  [7]  despite 
differences  in  the  kinetic  constants  used  between  the  two  studies.  Increased  U  increases  MMA 
concentration  as  Ts  increases.  CO2  concentration  increases  with  x(  as  more  of  the  product  is 

produced  along  the  plate  and  increases  with  U  due  to  increased  production  rate  as  shown  in 
Figure  5.  H2O  concentration  follows  the  CO2  profile.  Unlike  MMA,  CO2  and  H2O,  the  oxygen 
profile  is  significantly  different  and  is  shown  in  Figure  6.  The  oxygen  concentration  drops 
sharply  near  the  flame  attachment  point  to  zero.  Therefore,  most  of  the  PMMA  surface  sees  only 
the  monomer  vapor,  products  and  nitrogen  along  its  surface.  This  has  important  implications. 
Kashiwagi  et  al.  [22]  showed  that  PMMA  pyrolysis  kinetics  is  altered  by  the  presence  of  oxygen 
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in  flameless  pyrolysis  experiments  under  a  radiant  panel.  The  higher  the  oxygenconcentration 
the  lower  the  value  of  Ea  for  the  pyrolysis  kinetics,  which  is  attributed  to  O2  diffusion  into  the 
melt.  Our  calculations  appear  to  suggest  that,  under  the  boundary  layer  combustion  conditions, 
PMMA  pyrolysis  occurs  without  the  influence  of  oxygen  at  the  surface. 

The  NS  solutions  for  local  regression  rates  rQ  are  obtained  by  dividing  the  pyrolysis  rate  jfi, 

which  is  given  by  equations  (6)  to  (10),  by  the  polymer  density  and  is  plotted  in  Figure  7  for 
different  air  velocities.  The  regression  rate  is  determined  by  the  combustion  kinetics  and  by  the 
heat  transfer  from  the  flame  to  the  surface  due  to  convection  and  conduction  as  described  by 
equations  (6)  to  (10).  As  one  might  expect,  the  NS  solutions  are  less  accurate  at  extremely  small 
values  of  x]5  as  they  become  comparable  to  the  cell  size  (0.2  mm)  used  in  the  computations.  The 

regression  rate  is  the  highest  near  the  flame  attachment  point  due  to  the  small  distance  between 
the  flame  and  surface.  The  combustion  kinetics  plays  a  significant  role  in  determining  the  peak 
regression  rate.  Furthermore,  as  U  is  increased,  Da  decreases  and  flame  attachment  point  moves 
downstream  .  Therefore,  the  position  of  the  peaks  in  regression  rate  also  move  downstream  of 
the  plate. 

As  X!  increases  beyond  the  peak  shown  in  Figure  7,  the  kinetic  effects  become  less  significant 
and  transport  of  heat  is  the  controlling  mechanism.  As  xt  increases,  8  increases,  the  heat 
feedback  qw  to  the  surface  and  rQ  decrease  from  their  peak  values.  Furthermore,  as  the  air 
velocity  is  increased,  convective  and  conductive  transport  to  the  surface  increase  and  rQ  increases 

significantly  as  shown  in  Figure  7.  In  Figure  7,  the  area  under  the  curve  represents  the  total  mass 
loss  per  unit  time  per  unit  width  of  the  plate.  The  area  or  mass  burned  increases  steadily  with  xt 

and  reaches  70%  at  3  cm  from  the  leading  edge.  Clearly,  the  first  few  centimeters  of  the  plate 
contribute  significantly  to  the  total  mass  loss  rate  and  are  influenced  by  different  phenomena  than 
at  large  distances  downstream. 

3.2  Comparisonof  NS  with  PMMA  Experiments  for  Flame  Temperature 

Figure  8  shows  the  steady-state  NS  solutions  for  the  temperature  profiles  across  the 
boundary  layer  flame  at  different  locations  along  the  PMMA  plate  at  U=60  cm/sec.  The 
temperature  increases  with  distance  from  the  surface  as  one  approaches  the  combustion  zone, 
where  most  of  the  energy  is  released,  and  decreases  as  one  reaches  the  bulk  of  air.  Peak 
temperatures  reach  near  1900K  in  the  combustion  zone.  Near  the  leading  edge,  the  heat  release 
rate  due  to  combustion  is  high  and  the  peak  temperature  is  high.  At  large  values  of  x]5  the  flame 

loses  heat  to  the  ambient  by  convection,  conduction  and  radiation  as  the  boundary  layers  get 
thick,  and  the  peak  temperature  decreases  with  xr  Figure  8  also  shows  that  the  width  of 

temperature  peaks  increases  as  the  flame  thickness  increases  with  due  to  reduced  reaction- 

diffusion  rates.  Figure  8  also  shows  experimental  data  for  temperature  with  distance  from  the 
surface  y  and  at  different  distances  from  the  leading  edge  x(.  The  temperature  data  were  obtained 

about  1  min  after  ignition.  Figure  8  suggests  that  the  temperature  profiles  reached  steady-state 
quickly,  within  1  min.  Clearly,  Figure  8  shows  excellent  agreement  between  the  predicted  peak 
temperatures  and  the  data.  At  large  distances  from  the  leading  edge,  the  NS  solutions  seem  to 
over-predict  the  temperatures  on  the  fuel  side  of  the  flame  in  the  region  between  the  combustion 
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zone  and  the  surface.  In  the  experiments,  there  is  significant  soot  formation,  as  indicated  by  the 
brightness  of  the  flame.  The  soot  concentration  is  expected  to  be  high  on  the  fuel  side  of  the 
flame  and  increase  with  distance  from  the  leading  edge  due  to  increased  MMA  concentration. 
The  soot  enhances  radiative  loss  from  the  flame  and,  perhaps,  explains  the  slight  discrepancy 
between  the  data  and  the  computations  at  large  values  of  xr 

Figures  9  to  1 1  show  temperature  profiles  at  increased  values  of  U.  As  U  increases,  the  flame- 
to-surface  distance  S  decreases  and  the  gradients  increase.  Figures  9  to  1 1  show  reasonably  good 
agreement  between  the  numerical  predictions  and  the  data  at  the  increased  values  of  U.  At  large 
values  ofU=169  cm/sec,  the  model  appears  slightly  (about  1  mm)  to  under-predict  the  position 
of  the  temperature  peak  and  over-predict  the  peak  value  at  large  values  of  xr  In  the  experiments, 

the  flame  was  observed  to  flicker  slightly  at  large  x]  and  may  lead  to  lower  peak  temperatures. 

Furthermore,  y=0  represents  the  surface  of  unbumt  sample.  In  the  experiments,  the  surface, 
however,  regresses  and  the  measured  temperatures  are  much  higher  than  the  predictions  at  y=0, 
especially  at  small  values  of  x,,  as  one  might  expect. 

The  flame  is  located  at  the  peaks  of  the  temperature  profiles  shown  in  Figures  8  to  1 1 .  In 
boundary  layer  flames,  the  flame-to-surface  distance  5  increases  with  and  decreases  with  U  as 
indicated  by  the  positions  of  the  peaks  in  Figures  8-11.  The  classical  boundary  layer  solutions 
for  the  non-dimensional  8/Xj  given  by  equation  (15),  NS  solutions  and  the  experimental  data  are 
plotted  together  in  Figure  12.  The  NS  solutions  are  given  by  equation  (1 1).  Clearly,  NS 
solutions  and  the  data  show  that  8/x,  varies  both  with  Re  and  air  velocity  U,  unlike  the  classical 

boundary  layer  solution,  which  contains  a  single  parameter  Re.  The  NS  solutions  consider  finite 
rate  combustion  kinetics  and  exhibit  a  maximum  at  small  values  of  Re  before  approaching  zero 
as  Re  approaches  zero.  The  curves  for  U=  60  and  84  cm/sec  also  exhibit  peaks,  which  are 
outside  the  scale  shown  in  Figure  12.  As  expected,  the  classical  BL  solutions  approach  infinity 
as  Re  approaches  zero  due  to  flame  sheet  assumptions.  Furthermore,  at  a  fixed  value  of  Re,  8/xi 
decreases  with  increased  U  for  small  values  of  Re1/2.  As  U  increases,  the  “flame  attachment 
point”  moves  downstream  and  xi  increases  at  the  attachment  point.  Therefore,  the  peak  value  of 
8/xi  at  the  attachment  point  decreases,  and  8/xi  varies  both  with  Re  and  U.  As  xi  and  Re 
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Re  increase,  the  NS  solutions  run  parallel  to  the  classical  solutions  and  8/xi  appear  to  depend 
only  on  Re.  Furthermore,  the  NS  solutions  for  8/xi  remain  lower  than  the  classical  solutions.  As 
X|  and  Re  increase  further  approaching  the  post  flame  plate,  the  NS  solutions  deviate  from  the 
classical  solutions  due  to  the  lack  of  transpiration  cooling  on  the  inert  plate.  More  significantly, 
the  NS  solutions  for  8/xi  are  in  good  agreement  with  the  PMMA  data  for  different  values  of  Re 
and  U. 


3.3  Comparison  of  NS  with  Emmons  Theory  for  Burning  Rates 

Figures  13  to  16  show  the  dimensionless  regression  rate  Nu  predicted  by  the  NS  solutions 
given  by  equation  (1 1)  at  different  values  of  U  along  with  the  classical  theory  of  Emmons  given 
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by  equation  (14).  Both  the  NS  solutions  and  the  classical  solutions  show  that  Nu  increases  with 
Re  as  the  convective  transport  of  heat  from  the  flame  to  the  surface  increases.  However,  at  low 
air  velocities,  U=60  and  84  cm/sec,  the  increase  in  the  value  of  Nu  predicted  by  the  NS  solutions 
do  not  follow  the  square-root  law  unlike  Emmons  solutions  for  Re1/2  <60.  As  the  air  velocity  U 
is  increased,  the  NS  solutions  appear  to  follow  a  square-root  law  as  indicated  by  the  dark,  solid, 
relatively  straight  lines  in  Figures  15  and  16  even  at  relatively  low  values  of  Re1/2  >30. 
Therefore,  Figures  13  to  16  show  that  Nu  predicted  by  NS  solutions  contain  two  independent 
parameters  Re  and  U,  unlike  the  classical  solutions.  As  U  and  Re  are  increased,  the  axial 
diffusion  terms  become  smaller  and  the  NS  solutions  approach  the  classical  boundary  layer 
solutions  qualitatively.  Quantitatively,  however,  NS  predictions  are  higher  than  Emmons 
solutions.  For  example  at  U=168  cm/sec,  Re1/2=  80,  NS  solutions  predict  Nu=l  1  vs  6.4  (~40%) 
as  given  by  the  classical  boundary  layer  theory.  Indeed  for  U>120  cm/sec,  the  NS  solutions  may 
be  represented  by 


Nu  -  £  +  0.1  Re 1/2  (16) 


for  Re1/2>30.  The  slope  in  equation  (16)  is  found  to  be  identical  to  that  in  equation  (14).  Here, 
the  non-zero  intercept  “s”  depends  on  the  burning  rate  near  the  flame  attachment  point,  which 
depends  on  the  competition  between  convection  and  reaction  rates  represented  by  Da.  Indeed, 
one  can  obtain  the  value  of  s  by  fitting  equation  (16)  to  Figures  14  to  16.  The  intercept  increases 
with  U;  8=1.5  @  U=84,  s=2.0  @  U=120.0,  8=2.8  @  U=168  cm/sec.  This  may  suggest  that  e  is 
inversely  proportional  to  Da. 


Clearly  Nu  predicted  by  NS  solutions  are  higher  than  classical  BL  theory  due  to  various 
approximations  discussed  above.  This  is  qualitatively  consistent  with  the  calculations  of  Chen 
and  Tien  [6],  who  showed  that  the  NS  solutions  very  close  to  the  leading  edge  (Re1/2<30)  are 
slightly  higher  than  Emmons  as  shown  in  their  Figure  1 1 . 

3.4  Comparison  of  NS  with  PMMA  Experiments  for  Burning  Rates 

Figures  13  to  16  also  show  experimental  measurements  for  the  regression  rates  at 
different  values  of  the  bum  times.  In  these  experiments,  the  entire  sample  is  burned  for  a 
prespecified  time  t,  which  was  varied  from  5  to  20  min.  Clearly,  the  measured  burning  rates  vary 
with  t  significantly  unlike  the  theories  which  assume  steady  state.  Since  the  temperature 
measurements  shown  in  Figure  8  suggest  that  the  transport  processes  in  the  gas  phase  are  at 


steady  state,  the  long  transient  exhibited  by  the  data  in  Figures  13  to  16  must  stem  from  the 
unsteady  phenomena  in  the  condensed  phase.  Unlike  the  gas  phase  transport,  which  occurs 
rapidly  on  a  short  time  scale,  transport  in  the  condensed  phase  is  slow  due  to  high  thermal  mass 
and  occurs  over  long  time  scales  exhibited  by  the  data  in  Figure  13-16.  There  are  two  kinds  of 
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transient  phenomena  that  have  opposite  effects  on  Nu;  (1)  moving  boundary  or  curvature  effect  at 
small  values  of  x)  or  Re;  and  (2)  in-depth  pyrolysis  effect  at  large  values  of  x(  or  Re. 

3. 5  Moving  Boundary  Effects 

At  small  values  of  x)  or  Re,  Figure  13  shows  that  the  measured  Nu  at  t=5  min  are  slightly 

below  the  NS  predictions.  As  time  progresses,  the  deviations  between  the  NS  solutions  and  the 
data  increase.  The  steady-state  NS  solutions  and  Emmons  theory  assume  that  the  solid  surface 
remains  flat  as  it  bums.  In  the  experiments,  however,  the  initially  flat  surface  becomes  curved  as 
it  regresses  and  is  shown  in  Figure  17.  Figure  17  shows  the  formation  of  a  valley  or  a  step  near 
the  flame  “attachment”  point  as  time  progresses  and  is  discussed  by  Ndubizu  et  al.  [16,17,33-35] 
in  detail.  In  the  boundary  layer  configuration,  the  surface  regression  is  non-uniform  due  to  the 
variation  in  heat  feedback  and  rQ,  which  was  shown  in  Figure  8.  Figure  17  shows  that  the  surface 

curvature  increases  as  one  approaches  the  valley  from  the  leading  edge  of  the  plate  and  then 
decreases  with  Xj  as  the  heat  feedback  decreases.  The  variation  in  the  curvature  can  result  in  the 

formation  of  a  relatively  stagnant  region  behind  the  step.  Due  to  the  increased  barrier  for 
diffusion  and  increased  surface  area,  the  concentration  of  monomer  vapor  can  be  higher  and  that 
of  oxygen  lower  in  the  valley  than  on  the  flat  portion  of  the  surface.  This  can  lead  to  increased 
flame  stand-off  distance,  which  was  observed  from  the  measurements  of  temperature  profile 
across  the  valley  made  in  the  experiments.  Due  to  the  high  molecular  weight,  horizontal 
geometry,  short  length  of  the  fuel  plate  and  thin  and  viscous  melt,  very  little  melt  flow  was 
noticed  across  the  entire  length  of  the  plate  in  the  experiments.  The  flame  standoff  distance 
increases  with  time  as  the  valley  deepens  and  widens  as  shown  by  the  surface  profiles  in  Figure 
17.  The  increased  size  of  the  valley  causes  the  local  regression  rate  rQ  at  the  valley  to  decrease 
with  time  as  discussed  in  detail  by  Ndubizu  et  al.  [16,17].  Therefore,  the  data  for  Nu  shown  in 
Figure  1 3  decrease  with  time  unlike  the  steady-state  theories  at  small  values  of  x,  or  Re.  At  large 

values  of  x(  or  Re,  the  surface  is  relatively  flat  as  shown  in  Figure  17,  and  the  moving  boundary 
effects  on  Nu  are  negligible. 

At  t=5  min  and  small  values  of  Re,  the  circles  shown  in  Figure  13  are  close  to  the  NS 
solutions,  since  the  deviation  from  the  flat  surface  assumption  are  small  as  shown  by  Figure  17. 
Figure  17  also  shows  that  the  shape  of  the  solid  surface  deviates  from  being  flat  to  an  increasing 
distance  from  the  leading  edge  x|;  as  time  progresses.  Therefore,  the  deviations  between  the 

measured  and  predicted  values  of  Nu  extend  to  increasing  distance  Xj  or  Re  with  time  in  Figure 

13.  The  moving  boundary  effects  extend  to  about  Re1/2=  14  (or  x  =0.5  cm)  at  t=5  min  and  Re1/2= 

28  (x  =2.0  cm)  at  t=20  min  as  indicated  by  the  deviation  between  the  data  (circles,  diamonds)  and 

NS  solutions  for  Nu.  Indeed,  Figure  13  shows  the  spread  rate  of  the  moving  boundary  effects 
with  time. 

As  the  air  velocity  is  increased,  the  heat  transfer  to  the  surface  increases  at  a  given  value 
of  X!  due  to  increased  flame  attachment  distance  and  decreased  flame  standoff  distance  in  both 

the  experiments  and  theory.  Therefore,  increased  air  velocity  leads  to  increased  deviation  from 
the  flat  surface  at  a  fixed  bum  time  due  to  the  increased  heat  transfer.  This  leads  to  increased 
deviation  between  the  data  and  predictions  of  NS  solutions  with  the  increased  U  as  shown  in 


23 


Figures  14-16  at  fixed  values  of  the  bum  time  and  xr  For  example,  at  x=2.0  cm,  Figure  13 

shows  little  deviation  between  data  and  NS  solutions  at  U=60  cm/sec  (Rel/2=  28).  On  the  other 
hand,  Figure  16  shows  significant  deviation  at  x  =2  cm  and  U=168  cm/sec  (Re1/2=  46).  Figures 

14-16  also  show  that  the  moving  boundary  or  the  curvature  effects  on  Nu  spread  to  increased 
values  of  Re  with  increased  U,  at  a  fixed  value  of  t.  For  example,  at  t=  20  min,  Figure  16  shows 
that  the  moving  boundary  effects  extend  to  Re1/2=  65  (x=4  cm)  at  U=168  cm/sec  compared  to 

Re1/2=  28  (x=2  cm)  at  U=60  cm/sec  shown  in  Figure  13. 

At  large  air  velocities,  168  cm/sec,  both  NS  solutions  and  experiments  show  that  the 
flame  attaches  ~  3  mm  downstream  of  the  leading  edge  of  the  plate.  Within  the  quench  distance, 
spontaneous  combustion  could  not  be  sustained.  After  the  valley  or  step  was  formed  in  the 
experiments,  it  was  found  that  the  flame  and  the  valley  spread  upstream  very  slowly  toward  the 
leading  edge  of  the  PMMA  plate,  and  the  valley  gets  wider  as  the  time  progresses  [31,34,35].  It 
appears  that  the  valley  stabilizes  the  flame  and  enables  it  to  spread  upstream  across  the  quench 
distance.  As  the  flame  spreads  upstream,  a  fixed  point  on  the  surface  is  exposed  to  decreasing 
heat  flux,  which  can  result  in  decreased  regression  rate  with  time.  However,  the  upstream 
spreading  is  very  slow  and  is  expected  to  have  a  relatively  small  effect  on  rQ  compared  to  the 

effects  of  increasing  curvature  with  time. 

3.6  In-Depth  Pyrolysis 

At  large  values  of  x(  or  Re,  Figures  13  to  16  show  that  the  measured  values  of  Nu 
increase  with  the  increased  bum  time,  and  approach  the  NS  solutions  at  large  values  of  the  bum 
time.  This  is  opposite  to  the  effects  of  moving  boundary,  which  occur  at  small  values  of  Xj  or  Re, 

on  Nu.  The  NS  solutions  assume  that  the  pyrolysis  reaction  occurs  at  the  surface  rather  than  in 
the  bulk  of  the  polymer  melt.  Since  the  surface  has  no  volume,  it  reaches  steady  state  instantly, 
and  all  of  the  heat  flux  arriving  at  the  surface  is  expended  in  unzipping  the  polymer  as  described 
by  equation  (5).  In  the  experiments,  however,  the  heat  is  conducted  into  the  bulk  of  the  solid, 
which  takes  a  significant  time  to  melt  and  heat  up.  As  the  thermal  boundary  layer  gets  thicker, 
pyrolysis  occurs  at  the  surface  and  in  the  bulk.  Therefore,  only  a  part  of  the  heat  flux  arriving  at 
the  surface  is  used  to  unzip  the  polymer  in  the  experiments  at  any  given  time.  Experiments 
[3,10,1 1],  which  were  performed  by  exposing  a  PMMA  plate  to  a  uniform  radiant  flux  without 
combustion,  showed  that  the  pyrolysis  rate  increases  with  time  before  reaching  a  steady-state. 

The  transient  period  increases  from  less  than  5  min  at  5  W/cm  flux  to  20  min  at  1.6  W/cm  as 
shown  in  Figure  3  of  Reference  [3],  We  calculate  the  surface  heat  flux  from  Figure  8  using 
equation  (9).  The  surface  heat  flux  was  found  to  vary  from  3.2  W/cm  at  x  =0.4  cm  to  0.6  W/cm 

at  x=6.5  cm  for  the  conditions  shown  in  Figure  14  in  the  present  study.  We  expect  that  the 
transient  time  for  heat  up  to  be  significantly  small  in  the  present  study  since  the  surface  was 
melted  and  prepyrolyzed  prior  to  ignition.  Therefore,  the  leading  edge  region  exhibits  unsteady 
pyrolysis  due  to  the  moving  boundary  effects  rather  than  due  to  transient  heat  up  in  the  present 
study.  As  x,  or  Re  increases,  Figure  1 3  shows  that  it  takes  a  longer  time  for  the  pyrolysis  rate  Nu 
to  reach  steady  state  due  to  the  in-depth  pyrolysis.  For  example,  at  Re  =  34  (x  =  3  cm),  Nu 

reached  its  steady-state  value  within  20  min  as  indicated  by  the  small  deviation  from  the  NS 
solution  at  U=60  cm/sec  in  Figure  13.  However,  at  Re1/2=  50  (x  =  6.5  cm),  Figure  13  shows 
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unsteady  state  in  the  data,  which  are  significantly  lower  than  the  NS  solution  within  the  same 
time  period. 

Figure  13  also  shows  that  the  steady-state,  surface-pyrolysis  approximation  is  valid  for 
the  middle  region  of  the  PMMA  plate.  For  example,  the  data  represented  by  the  circles  are  close 
to  the  NS  solutions  between  Re1/2  =  14  (x  =  0.5  cm)  and  Re1/2  =  28  (x  =  2  cm)  at  t=5  min.  Below 

x  <0.5  cm,  moving  boundary  effects  are  significant  and  above  x  >2  cm,  in-depth  pyrolysis  effects 

are  significant.  As  time  progresses,  both  the  moving  boundary  effects  and  the  steady-state 
pyrolysis  front  spread  downstream.  The  data  represented  by  the  red  diamonds  are  close  to  the  NS 
solutions  between  Re1/2  =  28  (x  =  2  cm)  and  Re1/2  =  34  (x  =  3  cm)  at  t=20  min.  Thus,  Figure  13 

shows  the  spread  rates  of  the  moving  boundary  and  steady-state  pyrolysis  fronts,  which  are  an 
order  of  magnitude  lower  than  the  traditional  co-current  flame  spread  measurements  made  at  the 
polymer  surface  [42].  The  spreading  of  both  the  fronts  is  controlled  by  the  transport  in  the 
condensed  phase  rather  than  the  gas  phase  transport. 

As  the  air  velocity  U  is  increased,  the  heat  flux  to  the  surface  qw  increases  due  to  the 
decreased  8  at  any  given  xf.  Therefore,  pyrolysis  reaches  steady  state  faster  than  at  low  velocities 
as  shown  in  Figures  13-16  at  a  fixed  value  of  x,.  At  U=60  cm/sec,  Figure  14  shows  that  steady- 
state  was  not  reached  at  x=6.0  cm  (Re1/2  =  48)  in  20  min  as  indicated  by  the  deviation  between 

the  red  diamonds  and  the  NS  solutions.  At  U=168  cm/sec,  however,  Figure  16  shows  that  steady 
state  was  reached  within  the  20  min  at  x=6.0  cm  (Re1/2  =  80).  Figures  13  to  16  also  show  that 

the  steady  state  pyrolysis  front  spreads  to  increased  distance  Xj  with  increased  U  like  the  spread 

of  the  moving  boundary  effects.  For  example,  at  U=168  cm/sec,  Figure  16  shows  that  the  steady- 
state,  surface-pyrolysis  approximation  is  valid  between  x,=4  cm  (Re1/2  =  65)  and  x,=9.5  cm  (Re1/2 

=  1 00)  at  t=20  min.  On  the  other  hand,  at  U=60  cm/sec,  the  steady  state  pyrolysis  approximation 
is  valid  between  x=2  and  3  cm  at  t=20  min. 

The  increase  in  Nu  at  a  fixed  value  of  in  our  experiments  is  expected  to  be  different  in 

magnitude  from  the  flame  less  pyrolysis  studies  [10,1 1]  and  from  the  co-current  flame  spread 
studies  due  to  differences  in  the  initial  conditions.  In  the  flame  spread  [12,13,15]  and  flameless 
pyrolysis  studies  [10,1 1],  the  sample  was  at  room  temperature  at  t=0  unlike  in  the  present  study, 
where  the  surface  was  prepyrolysed  under  a  radiant  panel  prior  to  the  formation  of  the  boundary 
layer  flame.  Therefore,  the  surface  layer  was  melted  and  reaches  steady  state  quickly  similar  to 
the  study  of  Krishnamurthy  and  Williams  [2], 

3. 7  Comparisonsof  NS  with  Literature  Data 

V ery  few  studies  reported  regression  rate  data  for  the  forced  convection  boundary  layer 
flame  on  a  PMMA  plate  in  the  literature.  The  data  were  often  derived  from  co-current  flame 
spread  experiments  in  which  bum  time  were  not  clearly  defined  due  to  differences  in  ignition 
procedure.  Therefore,  it  is  difficult  to  make  quantitative  comparisons  among  the  data  and  the 
theories.  However,  qualitative  trends  may  be  observed  by  making  comparisons  between  the 
current  study  and  the  previous  data. 
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Agrawal  [15]  performed  experiments  on  co-current  flame  spread  underneath  a  horizontal 
PMMA  plate,  which  was  fixed  to  the  ceiling  of  a  rectangular  channel.  Air  was  forced  through 
the  channel  at  a  prespecified  velocity.  The  plate  was  ignited  by  a  pilot  flame  (3  cm  long)  near  the 
leading  edge  and  was  allowed  to  spread  to  the  end  of  a  76  cm  long  PMMA  (white)  plate. 

Agrawal  [15]  measured  the  surface  regression  at  different  values  of  x]  at  the  end  of  the  flame 

spread  experiment.  Agrawal  [15]  calculated  the  values  of  the  surface  heat  flux  from  the 
measured  values  of  regression  and  the  time  using  the  theory  given  by  Agrawal  and  Atreya  [3]  for 
the  unsteady  conduction  heat  transfer  into  the  solid.  Agrawal  and  Atreya  [3]  derived  analytical 
expression  for  regression  rate  as  a  function  of  time  for  specified  values  of  the  surface  heat  flux. 
The  surface  heat  flux  was  assumed  to  be  independent  of  time  since  the  gas  phase  transport 
reaches  steady  state  very  quickly.  Agrawal  [15]  calculated  the  values  of  the  surface  heat  flux  by 
substituting  the  measured  values  of  regression  and  the  time  into  the  analytical  expression  derived 
by  Agrawal  and  Atreya  [3].  Thus,  Agrawal  [15]  reported  the  steady  state  values  of  the  surface 
flux  profile  along  the  PMMA  surface  for  distances  between  0  and  40  cm.  The  surface  flux 
profiles  exhibited  peaks  with  Xj  and  were  found  to  be  significantly  different  from  Emmons 

solutions.  At  large  values  of  xp  however,  the  surface  flux  profiles  shown  in  their  Figure  2  in 
reference  [15]  seem  to  become  roughly  parallel  to  the  line  representing  a  square-root  relationship 
predicted  by  Emmons. 

The  surface  heat  flux  distribution  reported  by  Agrawal  [15]  are  non-dimensionalized 
using  equation  (1 1)  to  obtain  values  of  Nu  and  are  plotted  in  Figure  16.  The  NS  solutions  do  not 
extend  to  large  distances  unlike  Agrawal’s  [15]  experiments.  However,  Figure  16  appears  to 
show  reasonably  good  agreement  among  the  NS  solutions  (when  extrapolated),  our  experimental 
data  for  PMMA  for  large  bum  time,  and  Agrawal’s  [15]  data  for  80<Re1/2  <  140  (Xj  >  6  cm). 

Both  the  NS  solutions  and  the  data  appear  to  show  square-root  relationship  between  Nu  and  Re 
for  x  >6  cm.  For  x  <6  cm  (Re1/2  <  80),  Agrawal’s  [15]  data  are  well  below  the  NS  predictions 

like  the  data  obtained  by  Ndubizu  et  al.  [17]  at  large  bum  time  due  to  the  moving  boundary 
effects.  In  Agrawal’s  experiments,  the  first  few  centimeters  of  the  plate  might  have  developed  a 
large  valley  quickly  due  to  the  large  over-hang  (3  cm)  of  the  pilot  flame  during  ignition,  and  led 
to  large  moving  boundary  effect  on  Nu. 

Zhou  and  Femandez-Pello  [13]  conducted  experiments  on  co-current  flame  spread  over  a 
long  PMMA  plate  fixed  to  the  floor  of  a  rectangular  channel.  Zhou  and  Femandez-Pello  [13] 
reported  values  of  Nu  for  15  cm<x,<70  cm  and  performed  the  experiments  for  a  period  of  20 

min.  Therefore,  any  given  point  on  the  surface  burned  for  20  min  or  less  during  the  flame  spread 
process.  Since  no  measurements  were  reported  for  the  first  15  cm  of  the  plate,  one  might  expect 
little  effect  due  to  the  moving  boundary.  Figure  16  compares  the  NS  solutions  with  a  correlation 
of  the  data  reported  by  Zhou  and  Femandez-Pello  [13]  for  Re172  >125.  We  rescaled  the 
correlation  given  originally  by  Zhou  and  Femandez-Pello  [13]  to  Nu=0.076  Re172  based  on  the 
values  of  the  constants  A.w  =  0.052  W/mK  and  AT=1200  K  used  in  the  present  study  as  shown  in 
Figure  16.  Figure  16  shows  that  the  data  of  Zhou  and  Femandez-Pello  [13]  is  lower  than  the 
steady-state  NS  solutions  perhaps  due  to  significant  in-depth  heat  transport  at  large  values  of  xr 

Unlike  the  NS  solutions,  the  correlation  of  the  data  given  by  Zhou  and  Femandez-Pello  [13] 
seem  to  suggest  that  Nu  depends  only  on  Re.  However,  a  closer  look  at  their  data  given  in  Figure 
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5  of  reference  [13]  shows  a  slight  downward  deviation  from  the  correlation  for  Nu  at  large  values 
of  X!  and  could  be  due  to  large  effects  of  in-depth  pyrolysis. 

Knshnamurthy  and  Williams  [2]  formed  forced  convection  boundary  layer  flames  on  a 
PMMA  plate  (15  cm)  by  igniting  the  entire  surface  at  once  rather  than  by  flame  spread. 

However,  they  conducted  the  experiments  for  a  very  short  period  (~100  sec)  of  time. 
Knshnamurthy  and  Williams  [2]  showed  a  square-root  relationship  between  the  measured 
regression  rates  and  x,  at  high  oxygen  concentrations,  where  the  heat  flux  to  the  surface  might  be 

high  and  in-depth  pyrolysis  might  occur  relatively  quickly.  At  atmospheric  concentration  of 
oxygen,  however,  their  regression  rate  data  are  much  lower  than  Emmons  theory  and  are  similar 
to  the  data  presented  in  this  work  as  shown  in  Figure  16.  In  ambient  oxygen,  the  surface  heat 
flux  is  significantly  small  and  leads  to  large  transient  effects  due  to  in-depth  pyrolysis. 


4.0  CONCLUSIONS 

Solutions  of  NS  equations  show  that  Nu  depends  both  on  Re  and  U.  At  large  Re  and  U, 
the  solutions  approach  the  square-root  law  given  by  the  classical  Emmons  solutions.  Therefore, 
the  boundary  layer  approximations  significantly  affect  Nu  for  Re1/2<  120  considered  in  this  work. 
Unlike  the  theories,  PMMA  experiments  show  that  Nu  depend  on  Re,  U  and  bum  time  t. 
Comparisons  of  NS  solutions  with  the  experiments  show  that  the  steady-state  surface  pyrolysis 
approximation  is  valid  upstream  of  the  plate  at  short  bum  times  when  the  deviations  from  the  flat 
surface  are  small.  At  short  bum  times,  the  transient  effects  due  to  in-depth  pyrolysis  become 
increasingly  significant  at  increased  values  of  x,  and  the  NS  solutions  over-predict  the  data.  As 

the  bum  time  is  increased,  the  NS  solutions  over-predict  the  data  near  the  leading  edge  due  to 
increased  surface  curvature.  On  the  other  hand,  the  measured  values  of  Nu  increase  at  large 
values  of  xf  and  approach  NS  solutions  at  steady  state.  As  the  air  velocity  is  increased,  the 

transient  moving  boundary  effects  increase  and  the  effects  of  transient  in-depth  pyrolysis 
decrease  at  a  given  value  of  the  bum  time.  Therefore,  the  deviations  between  NS  solutions  and 
the  data  for  Nu  increase  upstream  and  decrease  downstream  as  U  is  increased.  When  the  time 
effects  are  taken  into  consideration,  both  the  NS  solutions  and  the  experiments  appear  to  be 
qualitatively  consistent  with  data  reported  by  Agrawal  [15],  Zhou  and  Femandez-Pello  [13]  and 
Krishnamurthy  and  Williams  [2],  Furthermore,  the  NS  solutions  for  the  temperature  profiles 
across  the  flame  are  in  good  agreement  with  the  data  for  the  entire  length  of  the  plate. 
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APPENDIX  A  -  CONSERVATION  EQUATIONS 
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Continuity 
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where  p  is  the  gas  density.  The  quantities  u  and  v  are  axial  and  vertical  velocity  components  in  x 
and  y  directions  respectively.  Here,  P  is  pressure,  E  is  total  energy  density  and  Ck  are  specie 
number  densities. 


A.l  Transport  Fluxes 

The  viscosity  p  effects  are  contained  in  shear  stress  x,  which  relates  to  velocity  gradients 
by  the  following  equations, 


A-2 
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The  conductive  fluxes  q  depend  on  gas  conductivity  A,  and  temperature  gradients  and  are  given 
by 


and 


dT  N 
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The  first  term  in  equations  (A9)  and  (A  10)  represent  thermal  conduction  ,  the  second  term 
represents  Defour  effects,  where  hk  is  specific  enthalpy  of  specie  k,  Av  is  Avagadro  number,  X  is 
gas  conductivity  and  T  is  gas  temperature.  The  diffusive  fluxes  in  x  and  y  directions  are  Udk  and 
Vdk,  respectively,  and  are  given  by 
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where,  Xk  and  Yk  are  mole  and  mass  fractions  of  specie  k  respectively,  and  Dkm  are  equivalent 
binary  diffusivities  of  specie  k  in  the  gas  mixture.  The  gas  mixture  contains  N  number  of 
components. 


A.  2  Thermodynamic  Equations  of  State 

Equations  of  state  further  define  the  total  energy  density  E,  specific  internal  energy  e  and 
specific  enthalpy  hk  as 


E  -  pe 


pe  = 
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and 


hk  =  a\k+a2kT  +  a3kT2  +a4kT\  k  =  l,2,...N,  (A15) 

where  alk...,a4k  are  polynomial  coefficients  for  specie  k  and  are  taken  from  sources  such  as  those 
compiled  by  NIST  (www.nist.gov),  University  of  Berkeley  (www.ucb.edu)  as  discussed  by 
Ananth  etal.  [28].  The  ideal  gas  law  is 

P=RTi ,  (A  16) 
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the  mass  specific  heats  Cpk  of  specie  k  are 


m,  $r’ 


k  =  1,2, ...N  , 


(A  17) 


and  the  mass  specific  heat  of  the  gas  mixture  Cpm  is 


=  ZYkCk 

*= i  k  pk 
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A.3  Transport  Properties 

The  transport  properties  are  calculated  from  kinetic  theory  of  gases  described  by 
Chapman-Enskog  equations.  The  diffusivities  of  specie  k  in  a  mixture  Dkm  are  calculated  from 
binary  diffusion  coefficients  D*j  and  mixture  composition  using 
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where  the  binary  diffusion  coefficients  are  related  to  temperature  by  fitting  a  power  law  in 
the  range  of  300  to  3000  K  to  evaluate  bfij  and  b2jj : 
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(A20) 


The  mixture  viscosity  is  calculated  from  specie  viscosities  and  mixture  composition  Yk  from 
the  following  equations 
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The  specie  viscosities  are  related  to  temperature  by  a  power  law  given  by 

ft  =cl,r2' ,  k  =  l,2,...N,  (A23) 

where  clk  and  c2k  are  known  constants.  Mixture  conductivities  A,m,  which  depend  on  the 
individual  conductivities  of  the  species  and  mixture  composition,  are  given  by 
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The  specie  conductivities  are  also  related  to  temperature  by  a  power  law 


Xk  =  d\kTd2‘  ,  k  =  1,2, ...N  , 


(A25) 
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where  dlk  and  d2|<,  are  known  constants.  Eucken  equation  is  used  for  the  conductivity  of  H2O 
vapor.  Therefore,  all  the  transport  properties  are  functions  of  both  the  temperature  and 
composition. 

A.  4  Radiative  Loss 

Heat  loss  Qr  from  the  hot  gases  to  the  ambient  is  assumed  to  follow  Stefan-Boltzman  law 

and  is 


q,=kpAt‘ -?:)■  (A26> 

The  gray  body  extinction  coefficient,  KPm,  is  calculated  based  on  the  concentrations  of  carbon 
dioxide  and  water  using 
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{All) 


where 

Kpk  =KX  +KJ  +  KJ2  +KJ 3  +K5T4  +K6T 5  +  K7T6  +  KJ1  , 

(A28) 


subscript  k  stands  for  CO2  or  H2O,  and  K)  to  are  known  constants. 

A.  5  The  BoundaryConditions 

The  above  partial  differential  equations  require  specification  of  boundary  conditions  on 
all  the  dependent  variables  p,  u,  v,  P,  E,  and  Ck  at  inlet,  outlet,  upper  and  lower  walls  of  the 
channel. 

A5.1  At  inlet,  x=0  and  aU  v  values ; 

The  inlet  velocity  components  uo,  total  energy  Eo, density  po  and  composition  as  number 
densities  Cko  are  specified.  These  boundary  conditions  assume  that  the  inlet  is  far  from  the 
combustion  zone  so  that  they  are  affected  by  the  changes  occurring  inside  the  computational 
domain.  The  inlet  pressure  P  is  unknown.  The  pressure  is  expected  to  change  along  the  channel 
due  to  momentum  boundary  layer  development  along  the  channel  walls  due  to  friction. 

However,  this  pressure  change  is  expected  to  be  small  compared  to  pressure  change  near  the 
leading  edge  of  the  porous  plate.  Therefore,  we  specify  that  the  inlet  pressure  drop  to  be  zero. 
They  are  written  as 

u(t,0,y)  =  u0,v(t,0,y)  =  0,E(t,0,y )  =  E0  ,  (A29) 
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A5.2  At  outlet  x=L  and  all  values  of_ y; 

Again,  we  assume  that  the  channel  outlet  is  far  from  the  combustion  zone  so  that  the 
outlet  conditions  are  not  affected  by  the  changes  taking  place  inside  the  channel.  We  set  all  the 
gradients  to  zero.  The  outlet  pressure  gradient,  however,  must  be  specified  to  obtain  unique 
solution  to  the  equations.  We  assume  the  outlet  pressure  to  be  atmospheric.  In  reality,  however, 
combustion  heat  is  released  all  the  way  to  the  end  of  the  channel,  affecting  the  outlet  conditions. 
But  the  heat  released  is  less  than  one  tenth  of  the  maximum  heat  release  rate.  These 
approximations  are  expected  to  be  accurate  for  high  fluid  velocities,  which  exist  in  forced 
convection  boundary  layers.  As  the  inlet  velocity  decreases,  they  become  less  accurate  and  the 
solutions  near  the  end  of  the  channel  (1-2  cm  near  exit)  are  not  accurate.  The  boundary 
conditions  are  written  as 


and 


du  _dv  _  dE  dp  dCk 
dx  dx  dx  dx  dx 


0 t ,  L,  y)  =  0, 


(A31) 


P{t,L,y)  =  Pam,  (A32) 


where  L=Li+L2+L3  . 


A5.3  On  top  solid  wall  and,  non-vorous  sections  1  and  3  of  the  bottom  wall  v=H,  v=0: 

The  top  wall  is  assumed  to  be  far  enough  from  the  combustion  zone  so  that  no  heat  is 
transferred  to  it  from  the  flame.  However,  there  is  no  slip  on  the  wall  and  no  mass  flow  through 
it.  They  are  written  as 


and 


dE 

dy 


dp  dCk 
dy  dy 


dP 

—  it,  x,  H)  =  0, 
dy 


(A3  3) 


u  -  v(t ,  x,H)  =  0 ,  (A34) 


for  all  values  of  x. 
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A5.4  At  impermeable  sections  of_ the  bottom  wall,  v-0 

The  bottom  wall  is  assumed  to  conduct  heat  away  from  the  flame  following  Newton’s 
law  of  cooling  with  a  characteristic  heat  transfer  coefficient  h.  As  a  result  of  this  boundary 
condition,  the  wall  will  heat  up  to  a  temperature  Tw,  which  varies  along  the  wall  and  is  calculated 
by  an  iterative  method.  The  adiabatic  and  constant  wall  temperature  boundary  conditions,  which 
were  used  in  the  existing  works  by  Mao  et  al.  [7],  Kodama  et  al.  [8]  and  Chen  and  Tien  [6]  occur 
as  special  cases  for  h=0  and  oo  respectively.  The  boundary  conditions  are  written  as 


and 


=  (A3  5) 

dy 


dp 

dy 


dC± 

dy 


^7"  0,  x,0)  =  0 , 
dy 


(A3  6) 


u  =  v(t,  x,0)  =  0 ,  (A3 7) 


for  x  <  L]  or  x  >  L1+L2. 


A.  6  Initial  conditions 

Initial  conditions  are  taken  as  the  inlet  conditions  at  x=0  over  the  entire  domain. 
Therefore, 


u( 0,  x,  y)  =  u0,  v(0,  x,  y)  =  0,  E{ 0,  x,  y)  =  E0,  (A3 8) 

and 

p{ o>  y)  =  Po  >  ck  (0,  ^  y)  =  cko  >  p(°,  x>  y)  =  pa,m  •  (A39) 

Initial  conditions  are  applied  prior  to  ignition. 
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